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We report a measurement of D^-D'^ mixing in D° — > K^-k^-k~ decays using a time-dependent 
Dalitz plot analysis. We first assume CP conservation and subsequently allow for CP violation. The 
results are based on 540 fb~^ of data accumulated with the Belle detector at the KEKB e"^e~ collider. 
Assuming negligible CP violation, we measure the mixing parameters x = (0.80 ± 0.29^o JJ'^ 1° 14)% 
and y — (0.33 ± 0.241q'5| lQ gg)%, where the errors are statistical, experimental systematic, and 
systematic due to the Dalitz decay model, respectively. Allowing for CP violation, we obtain the 
CPV parameters \q/p\ = 0.861°;^° ± 0.08 and arg(g/p) = (-14ti^l|t^-)o 

PACS numbers: 13.25.Ft, ll.30.Er, 12.15.Ff 



Mixing in the D'^-D^ system is predicted to be very 
small in the Standard Model (SM) jl| and, unlike in K°, 
and systems, has eluded experimental observa- 
tion. Recently, evidence for this phenomenon has been 
found in D" K+R- /tt+tt- and if+Tr" 

decays. It is important to measure D^-D'^ mixing in 
other decay modes and to search for CP-violating [CPV) 
effects in order to determine whether physics contribu- 
tions outside the SM are present. Here we study the 
self-conjugate decay D^^iCg tt+tt^. 

The time-dependent probability of flavor eigenstates 
and to mix to each other is governed by the 
lifetime Tjja = l/F, and the mixing parameters x = 
{mi — m2)/r and y = (Fi — r2)/2r. The parameters 
TOi,m2 (ri,r2) are the masses (decay widths) of the mass 
eigenstates |i:»i,2) = p\D^)±q\D^), and T = (ri-hr2)/2. 
The parameters p and q arc complex coeflticients sat- 



isfying \pY 



1. Various decay modes have 



been exploited to measure or constrain x and y [J]. For 



> Kg TT+TT decays, the time dependence of the Dalitz 



plot distribution allows one to measure x and y directly. 
This method was developed by CLEO using 9.0 fb 
of data; here we extend this method to a data sample 60 
times larger. 

The decay amplitude at time t of an initially produced 
\D°) or \D°) can be expressed as 



M{m^_,m\,t) 



A{rn?_, m^j^) 



ei{t)+e2{t) 



p ^2 
2 2 Nei(t) + e2(t) 



+ -A{m._,m^) , 

q -1-0 



(1) 



where A and A are the amplitudes for and \D'^) 

decays as functions of the invariant-masses-squared vari- 
ables m j_ = {Kg TT^ ) . The time dependence is con- 
tained in the terms ei^2{t) = cxp[— i(mj^ 2 ^ 2/2)^]- 
Upon squaring M and A4, one obtains decay rates con- 
taining terms exp(— Fi) cos(a;rt), exp(— Fi) sin(a;rt), and 
exph(l±2;)rt]. 

We parameterize the if^TT+TT" Dalitz distribution fol- 
lowing Ref. The overall amplitude as a function of 
771^ and m?_ is expressed as a sum of quasi-two-body am- 
plitudes (subscript r) and a constant non-resonant term 
(subscript NR): 

A{ml,ml) ^J2are'^'-Ar{ml,ml) + a^-^e'^^^, (2) 

r 

A{mi,m\) ^^CLrc'^- Ar{m\,ml.) + a^^^e'^^^ . (3) 



The functions A^. are products of Blatt-Weisskopf form 
factors and relativistic Breit-Wigner functions [7|. 

The data were recorded by the Belle detector at the 
KEKB asymmetric-energy e+e" collider 0. The Belle 
detector [9| includes a sihcon vertex detector (SVD), a 
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central drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), a barrel- like arrangement 
of time-of-flight scintillation counters (TOF) , and an elec- 
tromagnetic calorimeter. 

We reconstruct I?" candidates via the decay chain 
D*+ TT+D", £»" KItt+tt- 0. Here, tt^ denotes 
a low-momentum pion. the charge of which tags the fla- 
vor of the neutral D at production. The candidates 
are reconstructed in the 7r"'"7r~ final state; we require 
that the pion candidates form a common vertex sepa- 
rated from the interaction region and have an invariant 
mass within ±10 MeV/c^ of m„o- We reconstruct 

candidates by combining the candidate with two op- 
positely charged tracks assigned as pions. These tracks 
are required to have at least two SVD hits in both r-cf) and 
z coordinates. A Z?*+ candidate is reconstructed by com- 
bining the candidate with a low momentum charged 
track (the 7r+ candidate); the resulting £>*+ momentum 
in the e+e~ center-of-mass (CM) frame is required to be 
larger than 2.5 GeV/c in order to eliminate BB events 
and suppress combinatorial background. 

The charged pion tracks are refitted to originate from 
a common vertex, which represents the decay point of 
the . The D*^ vertex is taken to be the intersection 
of the momentum vector with the e'^e^ interaction 
region. The D'^ proper decay time is calculated from 
the projection of the vector joining the two vertices (L) 
onto the momentum vector: t = L ■ (p/p)(m^o/p)- The 
uncertainty in t (at) is calculated event-by-event, and we 
require crt < 1 ps (for selected events, {ut) ^ 0.2 ps). 

The signal and background yields are determined from 
a two-dimensional fit to the variables m„o and Q = 

(m^o —mj^o^^—m^)-c^. The variable Q is the kinetic 
energy released in the decay and equals only 5.9 MeV for 
D*'^ TrfD'-' decays. We parameterize the signal shape 
by a triple-Gaussian function for m^o^^, and the sum 
of a bifurcated Student t distribution and a Gaussian 
function for Q. The backgrounds are classified into two 
types: random tt^ background, in which a random tt^ 
is combined with a true D'^ decay, and combinatorial 
background. The shape of the m^o^^ distribution for 

the random tt^ background is fixed to be the same as 
that used for the signal. Other background distributions 
are obtained from Monte Carlo (MC) simulation. We 
perform a two-dimensional fit to the measured m^o^^-Q 
distributions in a wide range 1.81 GeV/c^ < m„o < 

1.92 GeV/c2 and < Q < 20 MeV. We define a smaller 
signal region \m^o^^ — m£)o| < 15 MeV/c^ and \Q — 

5.9 MeV| < 1.0 MeV, corresponding to 3cr intervals in 
these variables. In this region we find 534410±830 signal 
events and background fractions of 1% and 4% for the 
random tt^ and combinatorial backgrounds, respectively. 
The 771^0 and Q distributions are shown in Fig. [T] along 



with projections of the fit result. 
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FIG. 1: The distribution of (a) m^^o^^ with < Q < 20 MeV; 
(b) Q with 1.81 GeV/c^ < m^o^^ < 1.92 GeV/c^ Superim- 
posed on the data (points with error bars) are projections of 
the m,,a -Q fit. 

For the events selected in the signal region we perform 
an unbinned likelihood fit to the Dalitz plot variables 
?7i^ and rn\, and the decay time t. For decays, the 
likelihood function is 

1=1 j 

where j = {sig, rnd, cmb} denotes the signal or back- 
ground components, and the index i runs over D'^ candi- 
dates. The event weights f,- are functions of m„n and 
Q and are obtained from the m„o -Q fit mentioned 
above. 

The probability density function (PDF) 
T'sig (m?_ , , t) equals |7W(m?_, m+, t)p convolved 
with the detector response. Resolution effects in two- 
particle invariant masses are significant only for m^^. 
The latter, and variation of the efficiency across the 
Dalitz plot; are taken into account using the method 
described in Ref. The resolution in decay time t 

is accounted for by convolving V^^^ with a resolution 
function consisting of a sum of three Gaussians with a 
common mean and widths ak = Sk ■ iyt,i {k = 1 — 3). 
The scale factors Sk and the common mean are free 
parameters in the fit. 

The random tTs background contains real and 
D° decays; in this case the charge of the tTs is un- 
correlated with the fiavor of the neutral D. Thus the 
T'rnd PDF is taken to be (1 - ,f^)\M{m^_,m\,t)\^ + 
fw\M.{rn^_,rn\,t)\^ , convolved with the same resolution 
function as that used for the signal, where /w is the 
wrong-tag fraction. We measure /w = 0.452 ±0.005 from 
fitting events in the Q sideband 3 MeV< |(3-5.9 MeV| < 
14.1 MeV. 

For the combinatorial background, T'cmb is the prod- 
uct of Dalitz-plot and decay time PDFs. The latter is 
parameterized as the sum of a delta function and an 
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exponential function convolved with a Gaussian resolu- 
tion function. The timing and Dalitz PDF parameters 
are obtained from fitting events in the mass sideband 
30 MeV/c2< |m^o^^ - m^ol < 55 McV/c^. 

The hkelihood function for decays, £, has the same 
form as £, with M and M (appearing in T'sig and T'rnd) 
interchanged. To determine x and y, we maximize the 
sum In £-|-ln C. Table|I]lists the results from two separate 
fits. In the first fit we assume CP is conserved, i.e., 
a ~ a, (f) = (j), and p/q = 1. We fit all events in the 
signal region, where the free parameters are x, y, r^o, 
the timing resolution parameters of the signal, and the 
Dalitz plot resonance parameters ar(NR) and ^^(nr) • The 
fit gives r^o = (409.9 ± 1.0) fs, which is consistent with 
the world average The results for and (pr for the 
18 quasi- two-body resonances used (following the same 
model as in Ref. [6|]) and the NR contribution are listed 
in Table |TT1 The Dalitz plot and its projections, along 
with projections of the fit result, are shown in Fig.[2l We 
estimate the goodness-of-fit of the Dalitz plot through a 
two-dimensional test Q and obtain /ndf — 2.1 for 
3653—40 degrees of freedom {ndf) . We find that the main 
features of the Dalitz plot are well reproduced, with some 
significant but numerically small discrepancies at peaks 
and dips of the distribution in the very high m?_ region. 
The decay-time distribution for all events, and the ratio 
of decay-time distribution for events in the A'* (892)+ and 
if* (892)^ regions, are shown in Fig. [3] 

TABLE I: Fit results and 95% C.L. intervals for x and y, 
including systematic uncertainties. The errors are statisti- 
cal, experimental systematic, and decay-model systematic, 
respectively. For the CPV-allowed case, there is another so- 
lution as described in the text. 



Fit case Parameter Fit result 95% C.L. interval 

Wo ^?(%) 0.80 ±0.29in?io.M (0.0,1.6) 

CPV y{%) 0.33 ±0.24 tg-g^lg g^ (-0.34,0.96) 

CPV 0.81 ±0.30l° J?t° °^ |x| <1.6 

y{%) 0.37 ±0.251°:?^ tro8 1^1 <l-04 

\q/p\ 0.86l»i^iro'3±0.08 

arg(g/p)(°) -l^t^-l^.l 



TABLE II: Fit results for Dalitz plot parameters. The errors 
are statistical only. 



Resonance 


Amplitude 


Phase (°) 


Fit fraction 


if* (892)- 


1.629 ±0.006 


134.3 ±0.3 


0.6227 


A'o*(1430)- 


2.12 ±0.02 


-0.9 ±0.8 


0.0724 


7^2* (1430) - 


0.87 ±0.02 


-47.3 ± 1.2 


0.0133 


if* (1410)- 


0.65 ± 0.03 


111 ±4 


0.0048 


i<'*(1680)- 


0.60 ±0.25 


147 ± 29 


0.0002 


K*lm2)+ 


0.152 ±0.003 


-37.5 ±1.3 


0.0054 


7^0* (1430)+ 


0.541 ±0.019 


91.8 ± 2.1 


0.0047 


A"2(1430)+ 


0.276 ±0.013 


-106 ±3 


0.0013 


7!:*(1410)+ 


0.33 ±0.02 


-102 ±4 


0.0013 


A'*(1680)+ 


0.73 ±0.16 


103 ± 11 


0.0004 


p(770) 


1 (fixed) 


(fixed) 


0.2111 


w(782) 


0.0380 ± 0.0007 


115.1 ± 1.1 


0.0063 


/o(980) 


0.380 ± 0.004 


-147.1 ± 1.1 


0.0452 


/o(1370) 


1.46 ±0.05 


98.6 ±1.8 


0.0162 


/2(1270) 


1.43 ±0.02 


-13.6 ± 1.2 


0.0180 


p(1450) 


0.72 ±0.04 


41 ± 7 


0.0024 




1.39 ±0.02 


-146.6 ±0.9 


0.0914 




0.267 ±0.013 


-157 ±3 


0.0088 


NR 


2.36 ±0.07 


155 ±2 


0.0615 




(GeV^/c") (GeV^/c") 
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FIG. 2: Dalitz plot distribution and the projections for data 
(points with error bars) and the fit result (curve). Here, m^- 
corresponds to m^(_R'g7r+) for decays and to m^(ii'g7r+) 
for D'^ decays. 

CP-conservation case, and the values for the CPV pa- 
rameters are further constrained: jg/pj = 0.95+q'2o and 
arg((7/p) = (— 2+i[5)°. A check with independent fits to 
the and D " tagged samples gives consistent results for 
X (y): 0.58%±0.41% (0.45%±0.33%) and 1.04%±0.41% 
(0.21% ± 0.34%), respectively. 

We consider systematic uncertainties arising from both 
experimental sources and from the D'^ —*Kg Tr+Tr" decay 



For the second fit, we allow for CPV. This introduces 
the additional free parameters \p/q\, aTg(p/q), ar(NR) and 
0r(NR)- The fit gives two solutions: if {x, y, arg(p/g)} is 
a solution, then {— x, — y, arg(p/g)+7r} is an equally good 
solution. From the fit to data, we find that the Dalitz plot 
parameters are consistent for the and D° samples; 
hence we observe no evidence for direct CPV. Results 
for \p/q\ and arg(p/g), parameterizing CPV in mixing 
and interference between mixed and unmixed amplitudes, 
respectively, are also found to be consistent with CP con- 
servation. If we fit the data assuming no direct CPV, the 
values for x and y are essentially the same as those for the 
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FIG. 3: (a) The decay-time distribution for events in the 
Dalitz plot fit region for data (points with error bars), and 
the fit projection for the CP-conservation fit (curve). The 
hatched area represents the combinatorial background contri- 
bution, (b) Ratio of decay-time distributions for events in the 
A'*(892)+ and 74:*(892)" regions. 



model. We estimate these uncertainties by varying rel- 
evant parameters by their ztlcr errors and interpreting 
the change in x and y as the systematic uncertainty due 
to that source. The main sources of experimental uncer- 
tainty are the modeling of the background, the efficiency, 
and the event selection criteria. We vary the background 
normalization and timing parameters within their uncer- 
tainties, and we also set /w equal to its expected value 
of 0.5 or alternatively let it float. To investigate possi- 



ble correlations between the Dalitz plot (m 



-I-' ' 



.) dis- 
tribution and the t distribution of combinatorial back- 
ground, the Dalitz plot distribution is obtained for three 
bins of decay time; these PDFs are then used accord- 
ing to the reconstructed t of individual events. We also 
try a uniform efficiency function, and we apply a "best- 
candidate" selection to check the effect of the small frac- 
tion of multiple-candidate events. We add all variations 
in X and y in quadrature to obtain the overall experimen- 
tal systematic error. 

The systematic error due to our choice of — > 
decay model is evaluated as follows. We 
vary the masses and widths of the intermediate res- 
onances by their known uncertainties ll|, and we 
also try fits with Blatt-Weisskopf form factors set to 
unity and with no dependence in the Breit-Wigner 
widths. We perform a series of fits successively exclud- 
ing intermediate resonances that give small contributions 
(/9(1450), i^*(1680)+), and we also exclude the NR con- 
tribution. We account for uncertainty in modeling of 
the S'-wave tttt component by using K-matrix formal- 
ism [l2|. We include an uncertainty due to the effect of 
around 10-20% bias in the amplitudes for the K*{1A10)'^, 
/sTq (1430)+ and (1430)+ intermediate states, which we 



K^gTT+TT 



observe in MC studies. Adding all variations in quadra- 
ture gives the final results listed in Table [D 

We obtain a 95% C.L. contour in the {x,y) plane 
by finding the locus of points where —2\nC increases 
by 5.99 units with respect to the minimum value (i.e., 
— 2A ln£=5.99). All fit variables other than x and y are 
allowed to vary to obtain best-fit values at each point 
on the contour. To include systematic uncertainty, we 
rescale each point on the contour by a factor vl + r^, 
where is a weighted average of the ratios of systematic 
to statistical errors for x and y, where the weights de- 
pend on the position on the contour. Both the statistical- 
only and overall contours for both the CPF-allowed and 
the CP-conservation case are shown in Fig. ID We note 
that for the CPF-allowed case, the reflection of these 
contours through the origin (0,0) are also allowed re- 
gions. Projecting the overall contour onto the x, y axes 
gives the 95% C.L. intervals listed in Table HI After 
the systematics-rescaling procedure, the no-mixing point 
(0,0) has a value — 2Aln£ = 7.3; this corresponds to a 
C.L. of 2.6%. We have confirmed this value by generat- 
ing and fitting an ensemble of MC fast-simulated exper- 
iments. 



no Cf»V (Stat, only 
no CPV 

CPV (Stat, only) 
CPV 




-1 



X (%) 



FIG. 4: 95% C.L. contours for {x,y): dotted (solid) corre- 
sponds to statistical (statistical and systematic) contour for 
no CPV , and dash-dotted (dashed) corresponds to statisti- 
cal (statistical and systematic) contour for the CPl^-allowed 
case. The point is the best-fit result for no CPV . 



In summary, we have measured the D^-D^ mixing 
parameters x and y using a Dalitz plot analysis of 
Ks 7r+7r~ decays. Assuming negligible CP vi- 
and 



5+0.09+0.10\ 



olation, we measure x = (0.80 ± 0.29_q q7_q ;^4, 
y = (0.33 ± 0.24t°:°^t° °|)%, where the errors are sta- 
tistical, experimental systematic, and decay-model sys- 
tematic, respectively. Our results disfavor the no-mixing 
point X = y ^ with a significance of 2.2cr, while the 
one dimensional significance for a; > is 2.4(t. We have 
also searched for CPV; we see no evidence for this and 
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constrain the CPV parameters \q/p\ and sig{q/p). 
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